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ABSTRACT 

Exchange vector currents are calculated up to one-loop order (corresponding to 
next-to-next-to-leading order) in chiral perturbation theory. As an illustration of 
the power of the approach, we apply the formalism to the classic nuclear process 
n + p — > d + 7 at thermal energy. The exchange current correction comes out to 
be (4.5 db 0.3) % in amplitude giving a predicted cross section a = (334 ± 3) mb 
in excellent agreement with the experimental value (334.2 ± 0.5) mb. Together 
with the axial charge transitions computed previously, this result provides a 
strong support for the power of chiral Lagrangians in nuclear physics. As a by- 
product of our results, we suggest an open problem in the application of chiral 
Lagrangian approach to nuclear processes that has to do with giving a physical 
meaning to the short-range correlations that play an important role in nuclei. 
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1 Introduction 



Nuclear responses to electroweak probes are described by matrix elements of the cor- 
responding currents. The major contribution comes from one-body currents, specifically, the 
sum of the currents of each nucleon. The principal corrections come from exchange currents 
(or meson-exchange currents) that consist of two-body and higher-body currents. By now 
there exist a large number of unambiguous experimental evidence |lj for the presence and 
structure of meson-exchange currents. While phenomenological in character, the approaches 
taken so far can describe the large bulk of experiments || rather successfully, providing us 
our progressive understanding of nuclear processes. To the extent that nucleon-nucleon 
interactions are now fairly accurately understood, one can have a great deal of confidence 
in the theoretical tool with which the effect of exchange currents is calculated. There re- 
mains however the fundamental question as to how our phenomenological understanding of 
nuclear forces and associated meson currents can be linked to the fundamental theory of 
strong interactions, QCD. 

In the modern understanding of QCD, it is the spontaneous breaking of chiral symme- 
try associated with the light quarks that predominantly governs the structure of low-energy 
hadrons as well as the forces mediating between them. In fact, the full content of the 
gauge theory of strong interactions, QCD, can be expressed at low energy by a systematic 
chiral expansion starting with effective chiral Lagrangians ||. Stated more strongly, such 
an approach, known as chiral perturbation theory (ChPT), while reproducing the current 
algebra, is now considered to be exactly equivalent to QCD in long wavelength regime Q. 
The question that is immediately posed here is whether or not one can describe nuclear 
processes from the chiral perturbation theory point of view. 

Two recent developments suggest that the answer to this question is affirmative. The 
first is the work of Weinberg || and Ordonez, Ray and van Kolck [|j on understanding 
nuclear forces from chiral Lagrangians. The second is the explanation by the present au- 
thors J?] of the enhanced axial-charge transitions in heavy nuclei in terms of exchange axial 
currents in chiral perturbation theory treated to the same chiral order as for nuclear forces. 
Both of these two works provide a systematic and consistent field-theory-based understand- 
ing of the nuclear processes and validate at the same time the traditional approaches to the 
processes. In this paper, we present one more case where the chiral Lagrangian scores a 
stunning success, namely, the radiative neutron capture 

n +p — > d + 7 (1) 

at threshold. This result was briefly reported in a Letter ||]. Here we shall go into greater 
details of the calculation involved for the process ([!]). 

Historically, the process (||) was first explained quantitatively two decades ago by 
Riska and Brown ||] who showed that the ~ 10% discrepancy between the experimental 
cross-section and the theoretical impulse approximation prediction is eliminated by exchange 
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currents. Riska and Brown computed, using realistic hard-core wavef unctions, the two 
one-pion-exchange diagrams initially suggested in 1947 by Villars plus the uj and A 
resonance diagrams. That the dominant contributions to electroweak exchange currents 
could be gotten from current-algebra low-energy theorems was suggested by Chemtob and 
Rho 1 11] who gave a systematic rule for organizing the leading exchange-current diagrams 
effective at low energy and momentum. Although suspected since the Yukawa force was 
introduced, the work of Riska and Brown was the first unequivocal evidence for the role 
of mesons, in particular that of pions, in nuclear interactions. In this paper, we show that 
the terms considered by Riska and Brown are a main part of the terms that figure in chiral 
perturbation theory to next-to-next-to-leading (N 2 L) order and that when completed by 
the rest of the N 2 L order terms, chiral perturbation theory works impressively, confirming 
on the one hand the work of Riska and Brown || and on the other hand the earlier ( "chiral 
filter") conjecture of Kubodera, Delorme and Rho [jTsj| . 

In all cases studied so far, one is limited to long wavelength processes, with the typical 
energy /momentum scale Q much less than the chiral symmetry scale A x ~ 4-7T/V ~ 1 GeV. 
This is because ChPT is an expansion in Q/A x and its practical value lies where Q/A x -C 1. 
For nuclear physics, where both baryons and mesons enter on the same footing, the heavy- 
fermion formalism as discussed in []|, [l3j proves to be convenient for making a systematic 
expansion in derivative on pion fields as well as on baryons fields, d/A x and in m n /A x where 
m n is the pseudo-Goldstone boson mass. In this formalism, the nucleon mass is regarded 
as heavy, m N ~ A- l^ 1 



As in ||, the expansion is organized by the power v in Q" '. This expansion will be 
reviewed in the next section. Here we remark briefly how kinematical considerations affect 
the chiral counting rule, modifying it from the naive counting rule. The modification comes 
from the fact that the space part of the vector current and the time part of the axial-vector 
current (V 1 , A ) are proportional to (7*, 757 ) and hence are suppressed compared to the 
other part of the currents, (V°, A*). The correct counting for one-body currents in two- 
body systems is (V 1 , ^4°)ib = 0{Q~ 2 ). Compared to these one-body currents, the leading 
two-body current contributions are of order 0(Q) and the one- loop corrections are of order 
C(Q 3 ). Thus from the point of view of the chiral expansion, the two-body currents at 
one-loop order correspond to next-to-next-to-leading (N 2 L) order. This is the order that 
was computed in the case of the nuclear axial-charge transitions studied in and that 
will be adopted for the vector current matrix element for the process (|l]). Three-body (and 

#1 In principle, one can formulate chiral perturbation theory equally well in relativistic form for the baryons 
provided that care is exercised in arranging terms in consistency with chiral counting. At low nuclear matter 
density, however, the relativistic formulation is awkward and the heavy-fermion formalism is a lot more 
powerful and convenient in implementing in nuclear problems of the sort that we are concerned with. The 
situation will be different at higher matter density. This is because at higher matter density at which the 
nucleon mass is effectively reduced as discussed e.g., by Brown and Rho [O , H ], the heavy-baryon strategy 
may be suspect as the "1/mjv" expansion may not converge rapidly enough. One of the authors (MR) is 
grateful to Heiri Leutwyler for a discussion on this matter. 
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higher-body) currents are suppressed for exactly the same reason as for the suppression of 
3-body forces discussed by Weinberg ||]: the resulting 3-body currents are of order 0(Q 4 ) 
compared to the 1-body currents. 

Given these counting rules, the ChPT in nuclear systems can be made to correspond 
to computing Feynman diagrams involving external fields in the increasing power v embedded 
inside the most general process describing the transition from the initial nuclear state to 
the final nuclear state with interactions taking place before and after the current insertion. 
This strategy is essentially the same as what was first suggested, with somewhat ad hoc 
assumptions, in pl| . What this means in practice is that we are to take the most realistic 
nuclear wave functions and calculate the transition matrix elements with the ChPT graphs 
computed to the maximum possible order of chiral expansion. This allows us to fix the 
"counter terms" in the chiral expansion from experiments. At the present stage of our 
understanding on the role of chiral symmetry in nuclear systems, this seems to be the only 
practical attitude to take in confronting nuclear dynamics. We might point out that this 
point of view is consistent with the approach advocated by Weinberg || in his treatment of 
nuclear forces, in particular many-body forces. An important consequence of this strategy 
is then that only the current operators obtainable by ChPT are to be kept. This implies 
mainly two things. Firstly, it implies that short-wavelength effects are to be "filtered out". 
Secondly, it implies that ra-body currents with n > 2 are suppressed in nuclear systems 
with the mass number A > 2 in the same sense that n-body forces are suppressed ||. Of 
particular importance of the first implication is that when the currents are put in coordinate 
space, shorter-range interactions for ry± ^ r c ~ where r c is the hard core radius - cannot 
contribute in ChPT. This is because the "short-range correlations" that play an important 
role in nuclear observables represent those degrees of freedom that cannot be accessed by 
chiral perturbation expansion. Clearly this strategy would make sense only if the result does 
not sensitively depend on the precise value of the "hard-core" size r c , within the relevant 
range for application in nuclei to the order considered, say, A^ 1 < r c < (2m,,-) -1 . Note that 
this is roughly the range that can be reasonably described in ChPT for the NN potential 

1- 

This paper is organized as follows. In section 2, chiral perturbation theory as appro- 
priate for the problem at hand will be briefly summarized. The main purpose of this section 
is to define the notations and indicate the strategy specific to nuclear problems as outlined 
above and detailed in QJ. For a general and comprehensive review of chiral perturbation 
theory, we refer to recent review articles [jl6f| . Section 3 deals with the derivation of the 
exchange electromagnetic currents in nuclei to N 2 L order for vanishingly small momentum 
transfers. The formulas derived in Section 3 are applied in Section 4 to the classic nuclear 
process (|l]). Some further remarks and conclusion are given in Section 5. The Appendices 
give detailed formulas used in the main text. 
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2 Chiral Perturbation Theory 



Here we shall briefly review the specific aspect of chiral perturbation theory (ChPT) 
adopted for our calculation to one-loop order. This section is organized as follows. We start, 



in Section 2.1, with a rederivation of Weinberg's counting rules with a slight generalization. 
The relevant chiral Lagrangians in heavy-fermion formalism (HFF) will be written down in 
Section and the renormalization procedure within HFF in Section 2.3. The particular 



feature of ChPT in nuclear physics will be described in Section 2.4. 



2.1 Counting rules 

Here we rederive and generalize somewhat Weinberg's counting rule ||] along the 
line described in [f?J. Although we shall not consider explicitly the vector- meson degrees of 
freedom, we include them here in addition to pions and nucleons. Much of what we will 
obtain later without vector mesons turn out to be valid in their presence. In our work, we 
consider the vector-meson masses - which are comparable to the chiral scale A x - as heavy 
compared to the momentum probed Q - say, scale of external three momenta or m^. In 
dense medium, the situation is different as discussed in [|l^]. There the vector- meson masses 
can be considered small, in the sense of Georgi's vector limit. In our discussions, we do not 
encounter the dense regime in which vector-meson masses drop as a function of density. 

In establishing the counting rule, we make the following assumptions: Every interme- 
diate meson (whether heavy or light) carries a four-momentum of order of Q. In addition 
we assume that for any loop, the effective cut-off in the loop integration is of order of Q. We 
will specify what this means physically when we discuss specific processes, for this clarifies 
the limitation of the chiral expansion scheme. 

An arbitrary Feynman graph can be characterized by the number En(Eh,Ee) of 
external - both incoming and outgoing - nucleon (vector-meson, external field) lines, the 
number L of loops, the number 1^(1-^, Ih) of internal nucleon (pion, vector-meson) lines 
and the number of C disconnected pieces of the Feynman graphs. For an connected graph, 
for example, we should have C = 1. Each vertex can in turn be characterized by the 
number di of derivatives and/or of m n factors and the number (hi, ej) of nucleon (vector- 
meson, external field) lines attached to the vertex. Now for a nucleon intermediate state 
of momentum p 11 = m^v^ 1 + k^ where k^ = 0(Q) and is a constant 4- vector with 
~ (1, 0), we acquire a factor Q~ l since 

S F (m N v + k) = ^- = 0(Q~ 1 ), (2) 
v ■ k 

where we have used the propagator obtained by the HFF which will be explained in the 
next section. An internal pion line contributes a factor Q~~ 2 since 

A F (q 2 ;ml) = - T -L- T = 0(Q- 2 ) (3) 
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while a vector-meson intermediate state contributes Q° (~ 0(1)) as one can see from its 
propagator, noting that D F u (q 2 ;m 2 v ) = (-gT + ^f] A F (q 2 ;m 2 v ) = O(Q ) A F (q 2 ;m 2 v ), 



A F (q 2 ; m 2 v ) = T ~ T = O(Q ) (4) 

q A — m y —m v 

where my represents a generic mass of vector mesons. And a loop contributes a factor Q A 
because its effective cut-off is assumed to be of order of Q. Finally, the C separated pieces 
of the Feynman graph contributes a factor Q^ l ~ c ^ due to the (C — 1) times of the energy- 
momentum conservation factor, (2ir) 4: 5 4: Q2 i Pi) wnere (12iPi) denote the sum of incoming 
and outgoing momenta of a separate piece of the Feynman graph. 

We thus arrive at the counting rule that an arbitrary graph is characterized by the 
factor Q v with 

v = -I N - 2I V + 4L + 4 - AC + J2 d * ( 5 ) 

i 

where the sum over i runs over all vertices of the graph. Using the identities, 1^ + Ih + 
In = L + V — C where V denotes the total number of vertices of the Feynman graph, 
Ijj = ^(J2ihi — Eh), In = ^{J2i n i ~ En) and J2i e i = Ee, we can rewrite the counting 
rule 

We recover the counting rule derived by Weinberg || if we set Eh = hi = and E F = e.% = 0. 
With non-zero external heavy meson lines, it is generally hard to satisfy our condition that 
all the internal lines should carry momenta of order of Q. But there is no limitation on the 
number of hi, e% and E F - 

The quantity Vi is defined so that 

di + y + hi + e, - 2J > 0. (7) 

This is guaranteed by chiral symmetry || even in the presence of external fields. In the 
absence of external fields (as in nuclear forces), 

di + ^- + hi-2>0. (8) 

This means that the leading order effect comes from graphs with vertices satisfying 

di + ^ + hi -2 = 0. (9) 

Examples of vertices of this kind are: 7r k NN with k > 1 (di = 1, m = 2, hi = 0), 
hNN (di = 0, rtj = 2, hi = 1), four-Fermi contact interactions, (dj = 0, rjj = 4, hi = 0), 
hir k with k > 1 (di = 1, Hi = 0, /ij = 1), etc where h denotes vector-meson fields. In iViV 
scattering or in nuclear forces, = 2 and Eh = 0, and so we have v > 0. The leading 
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order contribution corresponds to v = 0, coming from three classes of diagrams; one-pion- 
exchange, one-vector-meson-exchange and four-Fermi contact graphs. In ttN scattering, 
= 1 and En = 0, we have v > 1 and the leading order comes from nucleon Born graphs, 
seagull graphs and one- vector-meson-exchange graphs. 

In the presence of external fields denoted generically £ , the condition becomes jn| 

dj + y+ /H-2j>-l. (10) 

The difference from the previous case comes from the fact that a derivative is replaced by a 
gauge field. The equality holds only for the vertices £tt j NN with j > (di = hi = 0, rii = 
2), Ehiti with j > (/ij = 1, Hi = di = 0) and £n k with A; > 1 (Tij = rtj = 0, di = 1). We 
will later show that this is related to the "chiral filter" phenomenon. The condition ( |l0|) 
plays an important role in determining exchange currents. Apart from the usual nucleon 
Born terms which are in the class of "reducible" graphs and hence do not enter into our 
consideration, we have two graphs that contribute in the leading order to the exchange 
current: the "seagull" graphs and "pion-pole" graphs both of which involve a vertex 



with V{ = 0. On the other hand, a vector-meson-exchange graph involves a Vi = 2 vertex. 
This is because di = 1, hi = 2 at the £hh vertex and di = 1, nj = 2, hi = 1 at the £hNN 
vertex, where £ denotes external fields. Therefore vector-exchange graphs are suppressed 
by power of Q 2 . This counting rule is the basis for establishing the chiral filtering even when 
vector mesons are present (see Appendix |C|). Thus the results we obtain without explicit 
vector mesons are valid more generally. 

We now focus on the nuclear responses to electroweak probes. Denoting the vector 
(V^) and axial-vector currents (A^) collectively by G^, we decompose the latter in terms 
of the number of nucleons involved: 

G" = G£ B + G& J + .». (11) 

In principle, there will be n-body currents for iV > n > 1 in iV-body systems. Compared 
to the leading 2-body currents, the leading 3-body currents are suppressed by a factor 
0(Q 2 /A 2 ,) from the counting rule eq.(||) and further suppressed by a factor 0(Q/A x ) for 
exactly the same reason as for the suppression of the 3-body forces |J. Thus the net 
suppression factor is of order (Q/A x ) 3 compared to the 2-body currents. We can therefore 
safely ignore three-body and higher-body currents. Limiting ourselves to the two-body 
subsystem, we have En = 0, S^- = 2 and Ee = 1. The counting rule reads 

v = l-2C + 2L + J2"i, Vi^di + ^ + hi + ei-2 (12) 

i 

where C = 2 for one-body currents and C = 1 for two-body currents. This counting rule 
would naively imply G^ B = 0(Q~ 3 ) and G^ B = 0(Q~ 1 ). But as mentioned in Introduction, 



#2 



These are standard jargons in the literature. See [Q, hl| 
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V i and A 


V° and A 1 


(naive) 


1-body 


0{Q- 2 ) 


0(Q~ 3 ) 


0(Q- 3 ) 


2-body (leading) 


0(Q~ l ) 


0(Q°) 


0(Q~ l ) 


2-body (1-loop correction) 


0(Q l ) 


0(Q l ) 


0(Q l ) 


(M 2B )/(M 1B ) 


OiQ 1 ) 


0(Q 3 ) 


0{Q 2 ) 



Table 1: The chiral counting for the vector and axial- vector currents for two-body subsys- 
tems. The last row gives the ratio of the amplitude of the two-body currents compared to 
that of the one-body currents. The last column shows the naive counting. 



there is a further suppression factor of order Q/A x for (V 1 , A°)\b and (V°, A l )2B, with 
i = 1, 2, 3. The suppression of one-body currents can be understood simply by recalling that 
the leading-order vector and axial-vector one-body currents are proportional, respectively, 
to u'y^u = (0(1), 0(Q/m N )) and ugAl^l^u = (0(Q/m N ), 0(1)) where u denotes the spinor 
of the nucleon. As for the two-body currents, the same suppression factor appears with the 
interchange of the role of the vector and axial- vector currents. The effective chiral counting 
for the currents is summarized in Table [I]. One can see immediately that the effect of 
the two-body currents can be sizable while the loop corrections are small for the space 
component of the vector current and the time component of the axial-vector current. This 



feature was noted sometime ago |12 



2.2 Effective chiral Lagrangian in heavy-baryon formalism 

The effective chiral Lagrangian that consists of pions and nucleons involving lowest 
derivative terms takes the form [18], 



1 A 



f 2 / \ f 2 

+ ^Tr (V M £tV^£j + ^m2Tr(S + £t) + • • • , (13) 



where 



A„ = liH^d^-^Hv. + A^ + y^-A^, 
V^S = ^E-i(V M + ^)E + fS(V M -^) (14) 

with the external gauge fields V M = V M • | and A^ = A^ • |. Here m N ~ 939 MeV is the 
nucleon mass, g A — 1-25 is the axial coupling constant, f n ~ 93 MeV is the pion decay con- 
stant and nix ~ 139 MeV is the pion mass. The ellipsis stands for higher derivative and/or 
symmetry-breaking terms that will be specified later as needed. Under chiral SU(2)xSU(2) 
transformation, the chiral field S = exp(i^) transforms as S — > g R 

SflJ G?r,Sl G SU(2)) 
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and the covariant derivative transforms the same as E does, 



V M E = a M E-i(V M + ^)E + iE(V M -^) 

-> gaVpEgl (15) 



with the external gauge fields transforming locally 



V M + ^ V^ + A^gniV^ + A^gi-id^-gl, 



The non-linear realization of chiral symmetry is expressed in terms of £ = v E = exp(i^j L ) 
and h = h(^, g h , g R ) defined with £ 

e - ffnC^ = Kgl 

The nucleon field A?" transforms as a matter field, i.e, N — > ZiiV, and its covariant derivative 



as D^iV -> /i-D^N while -> TiA^/i* wheref^ 



9,{J " ^'(V„ + AJ« - i«V„ - ^,){t. (16) 



Note that h(£, g h , g R ) is a complicated local function of £, g L and the explicit form of 
which is not needed. 

We have included the four-Fermi non-derivative contact term studied by Weinberg 
||]. We shall ignore possible four- Fermi contact terms involving quark mass terms and 
derivatives except for counter terms needed later. They are not relevant to the chiral order 
(in the sense defined precisely later) that we are working with. The explicit chiral symmetry 
breaking is included minimally in the form of the pion mass term. Higher order symmetry 
breaking terms do not play a role in our calculation. 

For completeness - and to define our notations - we sketch here the basic element 
of the heavy-fermion formalism (HFF) |H| applied to nuclear systems as developed by 



Jenkins and Manohar |2C] wherein the nucleon is treated as a heavy fermion. As stressed in 
Introduction, the relativistic formulation of ChPT works well when only mesons are involved 
but it does not work when baryons are involved since while space derivatives on baryons 
fields can be arranged to appear on the same footing as four-derivatives on pion fields, the 
time derivative on baryon fields picks up a term of order of the chiral symmetry breaking 



#3 We have defined two covariant derivatives involving chiral fields, V M S and A M . The first transforms as 
£ does while the latter as N does. We can express one in terms of the other, 

A M = j^(y^)^, V„E = 2£A„£ 

We use, as is done frequently in the literature, V M £ for the meson sector and A M for the meson-nucleon 
sector. 
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scale and hence cannot be used in the chiral counting. This problem is avoided in the HFF. 
To set up the HFF, the fermion momentum is written as 



p» = + p (17) 

where v^ 1 is the 4- velocity with v 2 = 1, and is the small residual momentum. (In the 
practical calculation that follows, we will choose the heavy-fermion rest frame = (1,0).) 
From the conventional nucleon field N(x), we define heavy fermion field B v (x) for a given 
four- velocity v^, by 

B v (x) = e Lm v v - x N(x), (18) 

which we divide into two parts, 

B v = flW + fit") = + ^-B v . (19) 

As defined, ^ +) (^ _) ) can be identified as positive (negative) energy solution. Since the 
free Lagrangian for the nucleon is 

N(i@ -m N )N = B v +) ifiB^+B^\i@ - 2m N )B v ~\ (20) 

B^ can be viewed as the relevant low-energy degree of freedom (with a zero effective 

>(-) 



mass), while B v as the irrelevant degree of freedom (with an effective mass of 2m N ). By 
integrating o 
fields. Since 



integrating out the Bv , we get an effective Lagrangian involving only the Bv and pion 



M +) =B< + \ (21) 

we can reduce all the 4x4 gamma matrices to 2 x 2 Pauli matrices. Define the spin operators 
Stf and S? by 

S£ = -IsM-f]. (22) 
S? = [S£,SZ]=ur a f } v a S vP . (23) 

In the rest frame, v» = (1, 0), so S° = = and 

\ ~ 2 ' 

gij = l e *iV, e 123 = l. (24) 

Denoting £?v by B and omitting the subscript v, our chiral Lagrangian ( |i~3|) expressed 
in terms of the heavy-fermion field to leading (i.e, zeroth) order in and in chiral order 
takes the form 

Co = B[ivD + 2ig A S-/\]B-\Y J C A (BT A Bf 

A 

f 2 f 2 
+ ^fTr ( V^StV^S ) + ^m*Tr(E + E t ) (25) 
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and the corresponding nucleon propagator SF(m N v + k) is 

S F ^ + k) = -^. (26) 

As mentioned, the HFF is based on simultaneous expansion in the chiral parameter and in 
"l/m N ". We have so far considered only the leading-order terms in l/m N , i.e., 0((l/m N )°). 
The -j— corrections can be gotten directly from the Lagrangian (13) pi], pj: 



= -^—B(-D 2 + (vD) 2 + S^F llu + g 2 A (vA) 2 + 2g A {vA,S-D})B (27) 
m N 2m N \ / 

where = f = - + [T^]. 

While eq.(p7|) is the first "l/m N " correction, it is not the entire story to the order 
considered. One can see that it is also the next order in the chiral counting in derivatives 
and expected in any case independently of the inverse baryon mass expansion. Generally, 
the coefficient appearing in each term should be viewed as a parameter rather than as a 
constant to be fixed by chiral symmetry in HFF. The most general chiral Lagrangian of 
order 0(Q) relative to the leading-order terms is 

A = ^BT 1/m B - A- (BB) (BiS ■ AB) - AL (BS»B) (SA, B ) (28) 

with 

T 1/m = -D 2 + {vD) 2 + {l + b l )S^T llv + {g 2 A + b 2 ){vA) 2 

+ 2{g A + b 3 ) {v ■ A, S ■ D} + 6 4 A • A + 2b 5 [v ■ A, S ■ D] 

+ b 6 m 2 n Tr(X + St) + b 7 ml (S + £*) _ l -b s S^(d^ - d u V%), (29) 

where denotes the external isoscalar vector field and fej (i = 1, 2, • • • , 10) denotes dimen- 
sionless constants. We have written down this part of Lagrangian for completeness although 
it turns out that none of the terms in ( |2"8| ) contribute to the process we are interested in to 
the chiral order considered. We have checked that there is no l/m N correction in our cal- 
culation. One might think that there would be non-zero contributions to V' 1 or A coming 
from the contact terms that have S ■ D or v ■ A, but these terms are forbidden by parity. 



2.3 Renormalization 

As in [fjffl, we shall adopt the dimensional regularization scheme to handle ultraviolet 
singularities in loop calculations. It has the advantage of avoiding power divergences like 
5(0) ~ A^ ui where A cu t is the cut-off mass. In d = 4 — 2e dimensions, all the infinities 
appearing in one-loop calculations are absorbed in i. To define our conventions and to 
apply the procedure to the case of the exchange vector currents, we repeat some of the 
essentials of ref . 0] . 
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Let us begin with the bare Lagrangian 

o o o o 

C = Co + Ci + £2 + 



(30) 



Here the subscript denotes the chiral power i/j as given byeq.(|). Note that the zeroth-order 

o I . o 

term Co is the same with the Lagrangian (|25|) with the replacement, B — > B = yZj^B, 

o 

7r a — ► 7r a = \fZ^-K a and all the couplings and masses replaced by their bare quantities. The 
wavefunction renormalization constants are 

3(M)9]AK) 



1 + 



4/ 2 



1 



2AK) 



3/ 2 • 

The renormalized quantities of interest are related to the bare ones by 



1 + 



u 

o 

9a 



A(ml] 



9A 



J n 
1 + 



dg 2 A 



A(ml 

w 

J - 



(31) 
(32) 



(33) 



where A(m 2 ) and 77 are, respectively, the quadratic and logarithmic divergences, 



A(mt 



m 



167T 2 
1 



4TTfJ, 2 



167T 2 



Te) 



mt 



4TTfl 2 



(34) 
(35) 



with [i the renormalization scale. 

For nuclear applications, we need to add four-Fermi interactions to (30). The relevant 
bare Lagrangian can be written as 



£4pt 

with 



o V I o \ 



o S 



Cs[BB\ +C s [BtB\ - ACt\BS^B\ -AC t [BS^tB 



oS 

C s 

oV 

C s 

oS 

Ct 

oV 



1 

-■x + 



1-d, d 2 -l 4 

"5,4 + , r 9 A 



1 - 3(d - 2) 5 2 
1 - 3(d - 2)4 



4 ™ 16 
A(m 2 



A(m 2 ) y 
J - 



f 2 
A(ml) 



f 2 

J 7T 



(36) 



(37) 
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We can easily see that the "counter terms" and Cj- are not necessary for removing di- 
vergences from one-loop corrections to the four- Fermi Lagrangian: There are no divergences 
involving S^. Such terms however could arise when one integrates out heavy-meson degrees 
of freedom. In the heavy-meson-exchange picture, the coefficients Cf and would come 
from the exchange of isoscalar mesons while Cg and C^f from that of isovector mesons. The 
coefficients C§ and C$ would come from mesons of even intrinsic parity (such as scalars 
and vectors) while Cj~ and Cj* would result from odd intrinsic parity mesons (such as pseu- 
doscalars and axial- vectors). In coordinate space, the Lagrangian (|36|) gives zero-ranged 
interactions, so in two-nucleon systems, only two terms effectively contribute. Weinberg|5|] 
chose Cg and (Cs and Ct in his notation) while in our calculations, we take C$ and 

There are many counter-terms that can enter in £2 at order 0(Q 2 ) compared to the 
leading chiral order terms. Some of them were written down in @ for the exchange axial 
currents. Here, we list only those terms that are needed for exchange vector currents: 

£2 = igA-p>~B S fl v u v a ( A^D u D a — D u A /I D a + D u D a A^) B 
+ 



7. 2 

A C7 Be^ VlI {A u ,T^}B 
+ i 9 -±czBv v v a SP[A v ,T aP \B 

+ {^^BAlV^T b va S a K ab B + h. C 

- GiBTf.uS^BBB-GiBT^BBS'^B 

- ^E^D a (Br a S^B)D^(Br a S^B) (39) 

where D^(BXB) = d^BXB) + B [X, rj B for any X, I> = TfY%, = d^F u - d u T^ + 
[T^,T U ] and 

1 2 i 

A a 6 = &ab - £ T a T b = ^ab ~ -Cabc^c, (40) 
V a(B = v a v p _ g a/3 _ ^gagfi ( 41 ) 

3 

All the divergences encountered in our calculation can be removed by 

d 3 o r? 

C2 = — — QaV + C 2 , 
C 7 = -ff + cH 



* 4 Due to the antisymmetry of the two-body wavefunction and since the force is zero-ranged, the contact 
Lagrangian can be rewritten as 

V* pt {r) = 6(r)[(C§-C$-C% -3C%) + 2P T (C% -C$)], P T = T1 ' ™ + 1 . (38) 

Thus, we can convert the Weinberg parameters (Cs, Ct) into our parameter (Cf , C$) by C§ = Cs — 2Ct 
and Cg = —Ct- 
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cs 

Gl + G2 

2 

G\ — G2 



R 



v 



2 V + c§, 

-A, 



l-(d-l)g\ + Afl(c s s -C v s +CS-CV)] +G§, 
l-4fi(c§ + C%-C$-C?) 



+ G% 



2/4 ' 



V.R 



(42) 



Note the C2 term has been introduced in M while the C7 and c% terms are new. They are 



consistent with Ecker's counter terms obtained by heat-kernel method [22]. The counter 
term eg is not associated with any divergence but is needed if the A-resonance degree of 
freedom is integrated out as we will do in our calculation. We also note that cf 1 , G^, R\ 
and eg are all finite, //-independent renormalized constants. 



2.4 Chiral perturbation theory in nuclear processes 

We wish to calculate operators effective in nuclei for transitions induced by the vector 
and axial- vector currents denoted respectively by and associated with the electroweak 
fields V M and A^. As mentioned, there are subtleties in applying ChPT to nuclear processes. 
Firstly, we need realistic initial and/or final nuclear wavefunctions. This is not a trivial 
matter in quantum field theory in general and in ChPT in particular. Perturbation series to 
bound states fail to converge because they consist of infinite number of "reducible" ladder 
diagrams which are of same chiral order. As stressed by Weinberg, chiral perturbation 
theory is useful in nuclear physics only for "irreducible" diagrams that are free of infrared 
divergences. This means that both in nuclear forces and in exchange currents, reducible 
graphs are to be taken care of by a Schrodinger equation or its relativistic generalization 
with irreducible graphs entering as potentials. This also implies that in calculating exchange 
currents in ChPT, we are to use the wave functions so generated to calculate matrix elements 
to obtain physical amplitudes. This is of course the standard practice in the theory of meson- 
exchange currents but it is also in this sense that ChPT is predictive in nuclei. We should 
stress that this precludes what one might call "fully consistent chiral perturbation theory" 
where nuclear forces, nuclear currents and wave functions are all calculated to the same 
order of chiral perturbation expansion. Such a calculation even if feasible is likely to make 
no sense unless one learns to solve QCD for nuclei. 

Let M be the transition matrix of the operator M 

M={V f \M\Vi) (43) 

where and \^f) are the initial and final wavefunctions, respectively. The A4 in ChPT 
is consistently expanded as 

M = M-2+ M-i+ M+i + M> (44) 
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where the subscript on M is v for Q u and M. > contains 0{Q 2 ) and higher in chiral ex- 
pansion. Our claim that will be substantiated below is that 7W> is small and hence can 
be ignored. This is guaranteed by the higher power in Q and by the "naturality" argu- 
ment. The power of ChPT is that we can actually calculate all first three terms accurately 
if we have accurate wave functions and accurate experimental information from other ex- 
periments, these to fix the counter terms appearing in the vertices. Since we ignore higher 
order terms, the terms retained fail to satisfy the exact continuity equation implied by the 
potential that gives rise to the accurate wave functions. But those "short-range" terms (and 
multi-body terms that figure in many-body systems) that must enter into W4> to satisfy 
the continuity equation is suppressed numerically by some order in the chiral counting and 
if the correction is small, then we are free to ignore them even though by so doing we will 
be violating certain Ward identities.^ If the chiral expansion reproduces QCD at some 
low energy and if the wave function is accurate enough, then the low-order expansion in 
the current must be arbitrarily accurate at that low energy. What we have accomplished 
in this paper is that with the wave function chosen, the chiral expansion converges for the 
thermal capture process, leaving errors of order 0.6 %! This means that the short-ranged 
term required by the continuity equation lies within that error bar. 

Let us consider many-body systems as in the case of axial-charge transitions To 
the next-to-next-to-leading order as in (|44]), chiral symmetry again says that at E <C A x , 
three-body and four-body etc currents and forces are zero. If one were to build an arbitrary, 
though phenomenologically guided, model unconstrained by chiral symmetry, there would 
be no reason why one should not include three-body, four-body etc. currents and forces 
and/or what might correspond to higher chiral order terms although he/she might find 
them small a posteriori. But then to satisfy one's continuity equation, one would have to 
keep adding a chain of what one might call "higher-order morass." What chiral symmetry 
says is that this is not an efficient way of doing physics if it has anything to do with QCD. 

In calculating loop graphs involving two-pion exchange, we encounter divergences 
that can be removed by four- Fermi counter terms and their second derivatives. When 
Fourier transformed into coordinate space, these terms appear as zero-ranged functions - 
<5 3 (r) - and their derivatives, all with renormalized coefficients. The prescription needed 
for handling these short-ranged terms does not follow from chiral symmetry alone and will 
have to be justified on a more general ground. As was explained in detail in Q, consistency 
with ChPT demands that those terms in exchange currents be (a) negligible in magnitude 
and (b) additionally suppressed by nuclear correlations when embedded in nuclear medium. 
To give an example, low-mass vector mesons (u, p) cannot contribute to the time part of 
the exchange axial current and the exchange magnetic moment operator while they do to 
nuclear forces. 



* 5 This means that we are satisfying Ward identities - associated with conserved currents - to the chiral 
order we retain and the violation will occur only at the next chiral order that is ignored. We note that this 
is a point frequently missed in the literature. 
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3 Exchange Currents 



Let us focus on the vector currents V ll (k;q) involving two nucleons with the kine- 
matics defined by 

ATi(pi) + N 2 ( P2 ) - Nx(p'i) + N 2 (p' 2 ) + V"(fc), 

where k^ is the momentum carried by outgoing external vector field and is the 
transferred four-momentum between nucleons (f" = \ (p 2 — p 2 + pi — p'\ ) . We are interested 
in the space component of the vector current with k^ — ► 0. The process (|l]) that we will 
analyze in detail involves the magnetic moment operator (MMO) 



1 

-V fc x J(fc;g) 



k=0 




where = V 3 '^ + -B^ is the baryonic (iso-scalar Lorentz-vector) currents. Hereafter, 
we will denote y 3 '^ simply by V*. We see from this formula that only terms linear in k^ 
can give nonvanishing contributions to the MMO. 

3.1 One-body (impulse) contribution 

In the long wavelength limit, the one-body vector current is dominated by the MMO 
which is entirely determined by experiments: 

e . . 

where \xs = fj, p + p» n ~ 0.87981 and /Uy = /x p — /x n ~ 4.70589. As indicated in Table 1, 
the corresponding space part of the one-body vector currents is of order 0(Q~ 2 ) in the 
chiral counting. The meaning of this operator in the context of ChPT is as follows. To 
leading order in the chiral expansion, the dominant contribution is the sum of single-particle 
transitions with vertices modified by chiral loop corrections. The renormalized constants 
(calculated to the chiral order we shall consider) can then be fixed by experiments. This 
means that there is no prediction involved at the level of the impulse approximation. Given 
the single-particle operator, prediction is made in two-body and higher-body corrections 
that can be computed in ChPT. 

3.2 Two-body contributions 

The Feynman graphs for two-body exchange vector currents are given in Fig. |. We 
shall now discuss their contributions in the ascending chiral order. 

3.2.1 Tree contributions 

It follows from the counting rule eq.(||) that the leading [0(Q~ l ) or 0(Q) relative to 
the impulse term) two-body currents come from the tree graphs, Fig. |2], with the vertices 
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Figure 1: Generic graphs contributing to the exchange vector current. The large filled 
circles represent one-nucleon and one-pion irreducible graphs and the solid (dashed) lines 
renormalized nucleon (pion) propagators. The graph (a) and (6) are the usual "seagull" 
and "pion pole" graphs with one-loop correction, respectively, and the graph (c) represents 
one-pion and one-nucleon irreducible graphs up to one-loop order. 



given by the bare ones (which will be renormalized at higher order and identified with 
the experimental values). Up to leading chiral order, the four-Fermi contact terms cannot 
contribute to the vector currents [o]. The leading two-body currents therefore are 



•( „ \z9a ^ 2 • <?2 . ,-, 
l{Tl X T 2 ) -f2 S l—2 J + (1 

J TT '"'TT 



4 



2), 



-2z(ti x r 2 ) 



, g\ S\ ■q 1 S 2 - q 2 



/2 ml -qf ml- q\ 



(47) 



where qi = p[ — Pi is the momentum transferred to the ith nucleon, q^ = —q^ — ^k^ and 
Q2 = ~ ■ The coupling constants and masses in the above equation are renormalized 
(measured) constants gA = 1.257 and f % = 93 MeV. The resulting MMO operator denoted 
Vtree, ( fr om eqs.(g7|I|)) 



far 



OA 171 * 



(48) 



/\/\/\!>- 



(«) 



(b) 



Figure 2: Tree graphs contributing to the exchange vector current, the "seagull" (a) and 
the "pion pole" (6) graph. 
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where 



and 



2/0 (x) 
Vi{x) 
V2{x) 

f(0) 
f(0) 



AtTX ' 

e~" 



Attx 
Attx 



(1 + x) 



-^yo(x), 



3 3 \ did 

1 + z + = x —z—yo( x > 



X X* 



dx x dx' 



(n 0t 2 ) z (<Ti 0<x 2 ), 



(Tl T 2 ) 2 



ff • (<Tl <X 2 ) - -(<Ti <X 2 ) 



(49) 



(50) 



with = ±, x. We should mention that this is the same as what one gets from the 
corresponding "pair" and "pionic" currents obtained in [pl|| . The only thing new here is its 
precise place in the chiral expansion. 



3.2.2 Loop corrections to the one-pion-exchange currents 

The first one- loop corrections corresponding to 0(Q 2 ) relative to the leading-order 
tree discussed above are renormalizations of the vertices in one-pion exchange graphs. We 
shall first show that the corrections appear only at the VttNN vertex. The ttNN proper 
function up to one-loop accuracy is M 



-IT,. 



a9A 



(51) 



where terms proportional to v ■ q are higher-order and hence can be neglected. There are 
thus no one-loop corrections to the ttNN vertex. Now the Ward identity restricts the Vtttt 
proper function to the form 



2ie 



3ab 



+ [q 



1 



2ie 3ab q» + 0(k 2 ) 



(52) 



where F^ 7rV (k 2 ) is the electromagnetic form factor of the pion, F^(k 2 ) = 1 + 0{k 2 ). This 
shows that there are no one-loop corrections from this vertex to the magnetic moment 
operator. Therefore we are left with the VttNN vertex appearing in Fig. [TJo,. Here the 
situation is quite analogous to the one-loop corrections in the exchange axial-charge currents 
P]. The one-loop Feynman graphs contributing to the relevant vertex are given in Fig. |3|. 
These one-loop graphs are computed in Appendix [A|. The result can be summarized simply 
as follows. The key point to underline is that for the magnetic moment operator, there 
are no corrections from the finite loop terms. The only contribution comes from the finite 
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(6) 



(/) 



/\/\/\^ .i 



I f 



(<0 



(5) 



Figure 3: One-nucleon, one-pion irreducible one-loop graphs contributing to the four-point 
V-kNN vertex. 



counter terms that describe the degrees of freedom that are integrated out from our low- 
energy effective Lagrangian. In fact there are two such terms: 



+ 



J 7T 



(53) 



where = ^fl 1 " and A M = ?f-A" Here we have defined the spin-| and isospin-| 
projection operators by 



T a Tb 



The resulting one-pion exchange current, denoted Vf^, is of the form 

V ln - lT 2 9f4 C 7 ™2 r ,2 € v vK a q^ 



9a Q2 ■ S 2 



4/4 C9 ml - q\ 

~ (§Tf + J (n x t 2 ) z ) (S? k a -g»k- S 1 )Vfq 2 p 
+ (l<->2). 



(54) 
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The counter terms bring in two parameters, and eg, to be fixed. These should be deter- 
mined from experiments. Relevant experiments would be isovector pion photoproduction 
or radiative weak amplitudes at low energy. However, the presently available data are not 
accurate enough to fix the constants sufficiently reliably. The difficulty lies in the fact that 
those processes are dominated by the seagull term and nucleon pole term and the constants 
we need to fix represent small corrections that cannot be extracted from the data. We shall 
therefore estimate them by the resonance saturation in a way analogous to what is done 



for the 0(Q ) Lagrangian counter terms for tttt scattering [23]. The relevant resonances are 



the uj and A saturating, respectively, the constants and eg. The results are 

- Cu) = = -0.1021, 

/I 8TT 2 g A (ml - ml) 

SA s g0J6(ff 
9f 2 9gA{mA - m N )m N 

where g^ is determined from the uj — ► n'y decay, g u = 17.55, and the iVA transition magnetic 
moment fiT and the irNA coupling C come from the fit to the A properties as explained in 
[24], C = —1.73 and = — 7.7 ± 0.5. For details, see Appendix 0. The resulting magnetic 



moment operator that we will use later is rather simple: 

»ln = tfn + f*h (56) 

with 



(57) 

This operator is equivalent to Fig. || calculated in the conventional meson-exchange current 
study. What is significant - and novel - in the context of chiral perturbation theory is that 
chiral symmetry constrains that there be no other contribution to the next-to-leading order 
one-pion exchange magnetic moment operator: There are no chiral loop corrections to the 
one-pion- exchange operator. This is in contrast to the case of axial-charge transitions (?J 
where chiral loop corrections to the one-pion exchange axial charge operator play a relatively 
important role. This suggests to lump all four terms of Fig. [2] and Fig. [I] together and call 
them "generalized tree operators". 

We should explain briefly why there are no genuine loop corrections ("chiral log 
terms") here in contrast to the axial-charge operator of J?]] and how the "counter terms" 
eq.(|53|) are saturated by the resonances. The genuine loop corrections arise from (in- 
dependent terms where q 11 is the momentum transferred by the pion. Now the nucleon 
propagators in a loop cannot bring in a q 2 -dependence because of the softness of the mo- 
mentum and v ■ q = 0. This means that if there is any momentum-dependent term, it must 
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A 



(a) (b) 

Figure 4: Graphs that correspond to the full one-loop corrections to the one-pion-exchange 
graphs. These graphs plus the tree graphs are referred to as "generalized tree" graphs. 

arise only from the momentum dependence of the pion propagators in the loop. We see from 
Fig. H that the first seven graphs (a — g) are momentum-independent and the next three 
graphs (h,i,j) are fe 2 -dependent but g 2 -independent. Therefore if any, chiral log terms can 
occur only at two-loop (or higher) order. 



3.2.3 Two-pion-exchange currents 

There are numerous diagrams of genuine loop character that can contribute in general 
kinematics to the two-pion-exchange two-body current. However a drastic simplification is 
obtained for the magnetic moment operator in heavy-fermion formalism: only four graphs 
(a), (b), (c) and (d) of Fig. || give non- vanishing contributions to the MMO. The graph (e) 
gives only zero-ranged operator and so will vanish in the matrix element, the graph (/) 
is non- zero but cannot contribute to the MMO as it does not contain terms linear in k^, 
the graph (g) is identically zero because it is proportional to v ■ S = and the remaining 
graphs ((h), (i), (j), (k) and (I)), being proportional to v^, do not contribute to the space 
component of the vector current. The calculation of the graphs is described in Appendix |B|. 

We pause briefly to explain why the zero-ranged graph (e) does not contribute in our 
calculation.^ Such zero-ranged two-pion-exchange graphs could be divergent (that is, they 
could have the pole). This divergence can be removed by a counter term generated 
by £2, leaving us a zero-range term with a finite coefficient. As in the case of the counter 
terms Cj' and eg, the (renormalized) coefficients should be fixed by experiments. However 
there are no experimental data, so we saturate them by resonances. As shown explicitly in 
Appendix [C], the coefficients obtained by this procedure come out to be simply zero. This is 
because the only low-lying resonance that can contribute is the p-meson, but its exchange 
does not give rise to terms linear in k^, and so gives vanishing contribution to the MMO. 
Even if the coefficients were non-zero, they would not contribute since short-correlations 
would "kill" them anyway unless the exchanged particle comes down below the chiral scale 
as in the case of the scalar meson in medium discussed in section 5. 

* 6 There are exceptional cases where this argument breaks down. This matter is discussed in section 5. 



20 




(0 (J) (*) (0 

Figure 5: One-pion and one-nucleon irreducible one- loop graphs contributing to the two- 
body vector current. Note that one-loop graphs which contain purely nucleonic intermediate 
state are excluded. Only the first four graphs (a, b, c and d) contribute to the magnetic 
moment operator. 
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The resulting two-pion exchange MMO takes the form 
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(58) 



The functions ko(x n ), ki(x w ) and hr(xn) are defined by 



ki(x n ) = 



k T (x n ) = 



1 



1 



r/ 3 



(2vr) 3 
d 3 q 



9 e i<?r 
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Analytic solutions are available for nonzero x n for the functions ko,i(x n ): 

1 



d_ 

dx 



k (x) = 
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h{x) = 



7TX^ 



irx' 



r Ki(2ar), 
(xK (2x) + 2Ki(2x) + xK 2 (2x)) 
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7TX 

— (K (2x) + 2xKi(2x)) 

7TX 



(59) 



(60) 



where K n (z) are the rath order modified Bessel functions. 

As we explained in several places, what we have to deal with is the finite-ranged 
contribution. For the n + p — > d + 7 process, the contributions proportional to the or 
Tj. + ^ are zero. Dropping zero-ranged operators that vanish by short-range correlations, the 
MMO can be expressed in terms of the "fundamental" constants qa and /„-, 
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(61) 
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where x n = m^r and the functions y's are defined in eq.(49). 



In summary, the total magnetic moment operator up to one-loop accuracy then is 

A* = MlB + Vtree + A*1tt + A^tt- (62) 

4 Application to n + _p — > d + 7 



4.1 A long-standing discrepancy 

The radiative np capture n+p — ► cZ + 7 at threshold is the most unambiguous nuclear 
process where the effect of pionic currents can be clearly seen. We shall apply the magnetic 
moment operator derived in chiral perturbation expansion to this process. 

The experimental cross-section measured accurately with the neutrons of thermal 
velocity, v n = 2200 m/sec, is 

o- exp = (334.2 ± 0.5) mb. (63) 

We first describe the prediction made in impulse approximation. For this, we note that the 
process involves very low energy, the relative momentum in the center of mass system being 

p = -m n v n ~ 3.4451 x 10 -3 MeV ~ (5.7278 x 10 4 fin)" 1 (64) 



2 

where the nucleon mass m N is defined by twice the reduced mass of the proton and neutron, 
m N = 2m p m n /{m p + m n ) ~ 938.92 MeV. At this extremely low energy, the process is pre- 
dominantly governed by the isovector Ml transition operator. We write the wavefunctions 
of the initial 1 Sq np state and the final deuteron state, respectively, by 

^ np (r) = 7= uo(r) XooGo, 
V47T r 



w(r) 



XUz Coo (65) 



'47r r 

where XJJ 2 (Ctt z ) is the spin (isospin) spinor and 

Si 2 (f) = 3 <ti • f <r 2 • r - cti ■ a 2 . (66) 

The radial functions uo(r) and u(r) are normalized so that their large-r asymptotic functions 
4>o(r) and cp(r) have (j>o(Q) = 0(0) = 1. We have 

hmnoW = Mr) = ^ 1 - ~, (67) 

r~*co sm(o s ) a s 



lim u(r) 

r— >oo 



hjnw(r) = ^(l + A + ^J (68) 
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where 77 is the D/S ratio, 7 = ^B d m N ~ 45.702 MeV ~ (4.3177 fin) -1 , B d ~ 2.2246 MeV 
is the deuteron binding energy, a s ~ —23.749 fm is the singlet np scattering length and 5 S 
is the 1 Sq np phase shift. In terms of the Ml transition matrix M defined by 

e 



2m p 



-M 



jd 3 r ttt p (r)[Ai(r)]*tt d (r;0), (69) 



the cross section can be written as 

a = GM 2 , G^ e Y A2s f s (70) 

where v n is the neutron velocity, u> — B d is the energy of the emitted photon, and A s ~ 
0.8846 fm" 1 / 2 is the normalization factor of the deuteron wavefunction, 

1 f f 00 r i 

= Jd 3 r *J(r)* d (r;0) = ^ * [n 2 (r) + w 2 (r)J . (71) 

Now the transition matrix Mib for the one-body magnetic moment operator ( |46| ) is 

POO 

M 1B = u v / (iruo(r)u(r). (72) 







Austern showed in 1953 [|26j that Mib can be estimated accurately in an effective-range 
expansion in terms of low energy parameters. The leading contribution is 



DC 



M = n v \ dr o (r)0(r) = ^ f 1 - — ) . (73) 



10 7 V 7«s 

The next order contribution in the effective range expansion can also be represented |26| by 
scattering lengths, effective ranges and the D-w&ve probability of deuteron. The impulse 
approximation obtained with this method comes out to be Gi mp (Austern) = (303 ± 4) mb. 
As we will explain in the following subsection, one can obtain a considerably more accurate 
value with the Argonne v is potential |£7| , 

Cimp = 305.6 mb. (74) 

The discrepancy between the impulse ( |74| ) (or rather the Austern result) and the exper- 
iment (^) had been a long-standing puzzle in nuclear physics until Riska and Brown || 
supplied in 1972 the missing 10 % contribution in terms of meson-exchange currents. In the 
following subsection, we will do an accurate one-loop chiral perturbation calculation, using 
the Argonne v±s potential. 

4.2 Exchange current contribution 

In calculating the matrix element of the two-body currents M2B, it is convenient 
to perform the angular integration first and evaluate the spin-isospin matrix elements for 
various spin-isospin operators: 



(^T^SS = /^(XOoClO^^fxioCoO 



(T^)sd = (xooCio) 1 TTt "^(fjxio Coo- (75) 



47T 



^(0) 
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It is straightforward to show that 

(T(- ) )ss = (Ti x) ) SS = 4, 



<4 -) >sd = = - 



16 



(76) 



and that all others are equal to zero. The resulting two-body transition matrix elements 
are 



M 2B = M 2 S | + M, 



rSD 
L 2B 



(77) 



with 
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where the functions fsT*\ x n) are ( x t = TO^r 
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4.3 The A degree of freedom 



So far, we have not considered the A degrees of freedom explicitly. It figured indirectly 
only as a resonance that saturates one of the counter terms. The role of the A is subtle 
because of the non-commutativity of the large- iV c limit and the chiral limit. The matter is 
simple if the mass gap, 5m = to a — m N , is either small enough to be regarded as 0(Q) or 
large enough to be regarded as 0(A X ). In the large- N c limit, 5m = 0(N~ 1 ) is to be taken 
to be smaller than m w = 0(N®). On the other hand, in the chiral limit, it is to be taken to 
be bigger than m v since it remains finite when m n goes to zero. In reality, neither is a good 
approximation. For the process we are considering, we believe that both approaches, one 
with the 5m treated as small (referred to as "type-I") and the other with the 5m treated 
as heavy (referred to as "type-II") could be employed equally well. 

In the type-I approach, we should treat the nucleon and the A on the same footing. 
Then the counting rule should be modified to 



v = A-2C- 



^ E N + E A 



+ E H + E E )+2L + J2^l v\ = d^ 



n,i + n* 



+ hi + ei -2 (80) 
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where is the number of external A lines and n* the number of A lines attached to the 
ith vertex. In this approach, we cannot integrate out the A degree of freedom without 
introducing non-locality in the effective Lagrangian. This also means that the A must be 
included in the "reducible graphs" or in terms of Schrodinger picture, in the wavefunctions. 
This procedure has often been used in nuclear physics calculations for processes that involve 
the excitation of the A resonance. Processes involving large energy transfer are better 
treated in this manner. 

In the type-II approach, we should be able to treat 5m as heavy, say, 5m = O(Q ). In 
this case, we cannot attach any Q-factor to the propagator, since l/(Sm — v ■ k) ~ 1/Sm = 
O(Q ). The counting rule in this case is 

u 11 = A-2C-(^ + E /1 + E H + E E \+2L+Y J vl\ »f = di+y +n*+/u+e;-2. (81) 

i 

If this applies, we can freely integrate out the A degree of freedom as well as other massive 
degrees of freedom. This will be the case for the process we are considering. As shown 
recently by Mallik |28|], this approach is equivalent to the type-I approach in certain kine- 
matics (such as in our case). It appears that the type-I approach is more general in the 
sense that it can be applied to a wide range of kinematical conditions whereas the type 
II is restricted. For our case, the type-II approach is much simpler. This can be seen as 
follows. Since the energy carried by the exchanged pion (since the photon is thermal) is 
small compared with the NA mass difference, the A propagator in tree graphs may be 
safely expanded as 

1 1 v ■ k , 

+ 7^9 + • • • . (82) 



5m — v ■ k 5m {5m) 2 
where v ■ k corresponds to the energy of the pion or the photon. In our one-loop calcu- 
lation, only the first term is retained. How good is this approximation? In general, the 
A propagator in the two-pion-exchange graphs may carry an arbitrary loop momentum. 
The only constraint we have is that the loop momentum be of order of the character- 
istic momentum, Q. In our process, however, the characteristic momentum involved is 



7 = yJBd m N ~ 45 MeV. Therefore the error involved in truncating the series ( p2| ) is at 
most of order r y/{5m) ~ 0.15 in the two-pion-exchange contributions. Since the overall 
two-pion-exchange contribution to the thermal neutron process is very small, as we shall 
find below, this is a safe approximation. 

4.4 Numerical results 

To calculate the capture cross section, we need accurate two-nucleon wavefunctions. 
This is because we need to have an accurate prediction for 1 Sq np scattering length, a s , to 
which the cross section is directly proportional. (There is also an implicit a s -dependence 
through the wavefunction.) A majority of nuclear potentials fail to reproduce the experi- 
mental value a e s xp = —23.749 ± 0.008 fm, because they do not fully incorporate the charge- 
independence-breaking (CIB) effect of the potential. It is not enough to incorporate the 
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r c (fm) 


0.1 
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Reid 


c(mb) 


336.1 


336.0 


335.7 


335.2 


334.4 


333.3 


331.8 


188.8 


Stree (%) 


2.608 


2.608 


2.609 


2.609 


2.605 


2.592 


2.563 


2.577 


st (%) 


0.435 


0.431 


0.423 


0.408 


0.382 


0.346 


0.301 


0.367 


st (%) 


1.364 


1.355 


1.334 


1.294 


1.228 


1.132 


1.011 


1.168 


61* (%) 


1.799 


1.786 


1.757 


1.702 


1.610 


1.477 


1.311 


1.535 


S2n (%) 


0.471 


0.457 


0.443 


0.424 


0.397 


0.362 


0.319 


0.319 


<5 2B (%) 


4.878 


4.851 


4.809 


4.745 


4.612 


4.431 


4.193 


4.431 



Table 2: The total cross section a and the 5's calculated with the Argonne v\% potential. 
The <5's are defined by S x = M x /Mib x 100 % for a given x, 5i n = 5f n + 5^ n , and 623 = 
<5tree + ^lvr + ^2vr- The last column shows the results with the Reid hard-core potential. 



Coulomb interactions and the mass difference between the charged and neutral pions. It is 
also necessary to incorporate the CIB effects in intermediate- and short-range potentials. 
Potentials which do not account for such CIB effects predict typically ~ a s ~ —17 fm. 
For example, Reid's hard core potential j2£§ predicts o s (Reid) = —16.7 fm, giving a poor 
impulse approximation result for the capture rate, <7i mp (Reid) = 173.7 mb to be compared 
with what is needed, ~ 306 mb (see (|83|)). A similar result is obtained with the Reid soft 
core potential, the Paris potential etc. An important point to note is that the ratios 5 of the 
matrix elements of the exchange currents to the impulse current are remarkably insensitive 
to the potentials, as we shall discuss later. 

We take the Argonne vis potential recently constructed by Wiringa, Stoks and Schi- 
avilla [27]. This potential is fit to 1787 pp and 2514 np scattering data in the range 0-350 
MeV with an excellent \ 2 of 1.09 and gives the deuteron properties - the asymptotic S-state 
normalization, As, the D/S ratio, rj, and the deuteron radius, dd - close to the experimen- 
tal values. Electromagnetic properties also come out well, modulo exchange-current and 
relativistic corrections. It predicts o* ft = —23.732 fm in excellent agreement with the ex- 
perimental value a e s xp = —23.749 ± 0.008 fm. The single-particle matrix element with this 



potential gives the impulse approximation cross section (given by fJ- tree of eq.(62)) 

Pimp = 305.6 mb, (83) 
rs -'9.6% less than the experimental value (Jexp — 

334.2 ± 0.5 mb. 

In computing the matrix elements of the two-body operators, we need to take into 
account short-range correlations in the wavefunctions. Short-range correlations must involve 
nucleon interactions at a length scale less than the scale given by chiral symmetry breaking 
and cannot be accounted for by low-order chiral perturbation expansions. In calculating 
Feynman diagrams of ChPT, the cut-off of the theory is incorporated in the regularization. 
Therefore nuclear wavefunctions do not incorporate the effect of the degrees of freedom that 
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Figure 6: Total capture cross section cr cap (top) and £'s (bottom) vs. the cut-off r c . The solid 
line represents the total contributions and the experimental values are given by the shaded band 
indicating the error bar. The dotted line gives Jtree, the dashed line <5trcc + 5^ K1 the dot-dashed line 
^trco + <5itt = <5trcc + 5^ + Sf^ and the solid line the total ratio, 623- 
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SD 


Sum (SS + SD) 


5 tree (%) 


1.534 
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1.228 
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0.449 
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0.397 


<5 2 B (%) 


1.864 


2.748 
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Table 3: The <5's with r c = 0.5 fm calculated with the Argonne v±s potential. 
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are cut off from the theory and the short-range correlation that is represented by a hard-core 
radius is missing. We do not know how to implement this effect in a way consistent with 
the strategy of ChPT. In this paper as in |7]], we account for the short-range effect simply 
by multiplying the integrand of the radial integral by 9{r — r c ), 



dr —* I dr . 

o Jr c 

This procedure is justified by the insensitivity of the result to the hard-core size. 

We use the physical values for masses and constants that appear in the theory. There 
are no unknown parameters other than the short-range cutoff parameter r c . The resulting 
total cross section is plotted and compared with the experiment [25] in Fig. |6] (top) for a 



wide range of r c , < r c <^ 0.7 fm. Figure || (bottom) shows the contribution of each term 
in terms of the ratios 5t re e = M tr ee/MiB, 6i n = M^/Mib, d~2iT = -^2?r/-^iB and 52B which 
is the sum. For 5i n , we divided it into two part, 5f w and df^. The ratio corresponding 
to the "generalized tree" contributions is Jgto = dtree + o~iw These ratios are plotted in 
Figure ||| (bottom). We also list the total cross section and the ratios in Table ||| and |||. 
In Table ||, the results are listed for varying r c 's. The last column shows the results with 
Reid's hard-core potential |29|] . Notice that, although the cross section is very sensitive to the 
potential, the ratios are essentially the same. This confirms that the ratios are highly model- 
independent, a conclusion reached previously.^ 7 We should also point out that the r c which 



gives approximately the same result with Reid's potential is r c ~ 0.6 fm to be compared 
with the hard-core radius of the Reid potential, 0.55 fm. In Table ^, we decompose the 
ratios into the SS and SD channels (where SS (SD) denotes the contribution to the S-w&ve 
(D-wave) of the deuteron wavefunction) for a fixed r c , r c = 0.5 fm. The intrinsic uncertainty 
associated with short-distance physics notwithstanding, the theoretical prediction 

o th = 334 ± 3 mb (84) 

(where the theoretical error bar represents only the dependence on the hard-core cut-off) is 
in remarkable agreement with the experiment, say, within less than 1% ! We have checked 
that the isoscalar moment contribution to the process is totally negligible. ^ 



4.5 The chiral-filter phenomenon 

One of the principal results in our calculation is that the "generalized tree" contri- 
butions dominate to the N 2 L order, with only a small correction (less than 0.6 % of the 

#7 What this means as far as the np capture is concerned is that all that is needed for an agreement with 
the experimental capture rate is to fix parameters of finite-order chiral perturbation theory so as to fit the 
singlet scattering length a 3 , since the ratios of the matrix elements are insensitive to nuclear potentials. 
It thus appears that such a refined potential as the Argonne potential is not really necessary for a 
satisfactory result. We are grateful to Jim Friar for emphasizing this point. 

* 8 A rough estimate using Figure 7(d) as the isoscalar exchange-current operator shows that the isoscalar 
contribution to the cross section is suppressed relative to the isovector one by a factor of more than 10 5 . 
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single-particle matrix element) coming from the genuine one-loop correction. This agrees 
with the "chiral filter" mechanism seen in the axial-charge transitions 0] and confirms the 
conjecture made in We should make a few additional remarks on this point. 

It is interesting to compare the Ml transition (thermal neutron capture) with the 
axial-charge transitions (7j. In fact, the meaning of "chiral filter" is slightly different in the 
two : In axial-charge transitions, the chiral filter mechanism is reflected in the smallness of 
the loop contribution compared with the tree contribution whereas in Ml processes, it is in 
the suppression of the two-pion exchange (one-loop) term compared with the total one-pion 
exchange (generalized tree) contribution. If we denote the operator by Ai A {Ai v ) and its 
matrix element by M A (M v ) for axial-charge (Ml) transitions, then the statement is 



M, 



loop 



ML + Ml « M 



M 
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2tt 



'Itt 
rV 



L treei 



« M v GTO = MY ree + M{„. 



(85) 
(86) 



These two observations can be unified if one accepts that what chiral filter conjecture really 
says is that the long-range operator (characterized by its tail e _m?rr ) dominates whenever 
it is not suppressed by kinematical or symmetry reasons. In the axial-charge transition 
case, as was discussed in the Aif^ can be decomposed into two parts, long-ranged and 
short-ranged as 



SsoftM A . ee + M l7rshort 



(87) 



with 

$soft 



C 3 f2 + 
J - 



7T \ , o , / 17 7T 

^)-0 + ^)( T -^ 



0.051 



where ■M.u z . sfiart represents the short-range contribution (shorter-ranged compared to one- 
pion-exchange) and the counter term cf is given by the charge radius of the proton. In (7j, 
we showed that 5 so ft (or c§) could be understood in terms of the p-meson degree of freedom. 
Since two-pion-exchange contributions are short-ranged, we can clearly divide the exchange 
axial-vector current into a long-ranged component and a residual short-ranged part: 

M short = Mi n . short + M 2n - 



Put in this form, the chiral filter effect is even more transparent and striking 

A 

~ (0.03 - 0.05) 



M s hort 
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ong 



(90) 



with the matter density ranging < p < po- I n the Ml transition in question, the corre- 
sponding long- and short-ranged operators are 



Ml ng 
J 1 short 



MY ree + MX^ + M\^ A = M GTO , 

Ml. 



(91) 



We see that in this extended form the chiral filter works equally well for both processes. 
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5 The BR Scaling 



There is very little understanding of the meaning of the short-distance correlation 
cut-off r c in the context of ChPT. As discussed in |7|, the loop terms contain zero-range op- 
erators in coordinate space. In addition, four-Fermi counter terms in the chiral Lagrangian 
with unknown constants are of contact interaction. At higher chiral order, increasingly 
shorter-ranged operators would enter together with the zero-ranged ones. Now if we were 
able to compute nuclear interactions to all orders in chiral perturbation theory, the delta 
functions in the current would be naturally regularized and would cause no problem. Such 
a calculation of course is an impossible feat. The practical application of chiral perturba- 
tion theory is limited to low orders in the chiral expansion, so a cut-off would be needed 
to screen the interactions shorter-ranged than accessible by the chiral expansion adopted. 
This means that all interactions of range shorter than the inverse chiral scale ~ 1 GeV -1 
should not contribute as additional operators. Clearly such a calculation would be mean- 
ingful only if the dependence on the cut-off were weak. Our calculation here meets that 
criterion. In the present work, we find that the r c -dependence is trivially small over wide 
range, < r c <^ 0.5 fm, and the theoretical prediction agrees exactly with the experiment 
at r c ~ 0.5 fm. The r c -independence and the success can be traced to the fact that there 
is no relevant low-energy resonances which contribute to the counter terms associated with 
four-Fermi interactions and the fact that the short-range correlation effect is automatically 
included in the np and deuteron wavefunctions. 

There appear subtleties, however, when nuclear density goes up, as in the case of 
the axial-charge transitions in heavy nuclei studied elsewhere. As nuclear density increases, 
some of the massive degrees of freedom that have been integrated out could come down 



in energy as suggested by Brown and Rho[14, 15]. If the mass comes down below the 



relevant chiral scale, roughly equivalent to r~ , then the short-range correlation can no 
longer suppress such a contribution. We shall now suggest that this is what happens with 
a scalar "dilaton" field. 

Consider the scalar field x associated with the trace anomaly of QCD exploited in 



[14, 15]. In matter- free space (p ~ 0), the scalar field, associated with scalar glueball, is 
massive with a mass m x ~ 2 GeV > m p , so may be integrated out. A simple calculation 
shows that this leads to an 0(Q 3 ) four-Fermi counter-term Lagrangian of the form 

AC = a(BB) (BY 1/m B) (92) 

where a is a dimension-4 constant and Tij m is the operator which appears in l/m N expan- 
sion, eq.(p9]), 



T 1/m = -D 2 + (v • Df + S^T^ + g 2 A (v • A) 2 + 2g A {v • A, S • D} + • • • . 

One can see from dimensional counting that this Lagrangian is highly "irrelevant" and 
hence would be negligible. However the situation can be drastically different in medium. 
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As discussed in 15], a low-lying excitation with the quantum number of scalar meson 
exists in dense medium. Such a "dilaton" is supposed to play an important role in nuclear 
physics as manifested in Walecka's theory with a low-lying scalar meson and also required 
by a "mended symmetry" emerging at a higher energy- momentum scale [pCj] . What this 
means is that instead of the zero-range interaction coming from the four- Fermi interaction, 
we would have a finite-range interaction with the range determined by the scalar mass in 
medium, m* where a is the scalar corresponding to the dilaton. We can gain a simple idea 
by noting that the Lagrangian ( |92|) can be obtained from the "l/m N " Lagrangian (|28| ) by 
replacing the nucleon mass by its effective operator, 

where po is the normal nuclear matter density and k is a constant. If one assumes that k is 



saturated by the a meson effective in medium [15], then one finds 



sIpo (94) 



m N ml 



where g a is the aNN coupling and m a the a mass. In mean-field approximation, (^) and 
( |94]) give, for m a ~ 600MeV and g a ~ 10, a substantial enhancement to the one-body axial 
charge operator as well as to the one-body isovector magnetic moment operator: 

5 a A° 5 a n p p 

7TT = = n — ~ 0.4 — . (95) 

A° fj, po po 

This is basically the enhancement in the axial charge in heavy nuclei as found in phenomeno- 



logical models with a low- lying scalar meson by Towner |3l]] and by Kirchbach et al [32]. 
It is equivalent to the Brown-Rho scaling f*/f-w ~ rn* N /tun as discussed in j^] . 



6 Conclusion 

This paper provides a chiral symmetry interpretation of the long-standing exchange- 
current result of Riska and Brown. It supports the chiral-filter argument verified in the axial- 
charge transitions and offers yet another evidence that chiral Lagrangians figure importantly 
in nuclear physics. 

Here as well as in |7|, ChPT is developed in the long wavelength regime where the 
expansion makes sense. The question remains as to at what point the expansion of ChPT 
breaks down. In particular, the meaning of short-range correlation remains to be clarified, 
an issue which cannot be ignored in any nuclear process. An intriguing question in this 
connection is posed by the observation that the generalized one-pion exchange process 
dominates even in e + d — > n + p + e' at a large momentum transfer |j|]. The only way 
one can understand this phenomenon would be that the chiral filter takes place even to 
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large momentum transfers. It would be interesting to verify this in a systematic chiral 
perturbation calculation. 

Another intriguing problem is the role of many-body interactions from chiral symme- 
try point of view. We know that for processes involving small energy transfers, three-body 
and other higher-body currents are suppressed. But this will not be the case for non- 
negligible energy transfers. There multi-body currents will certainly become important. 
Furthermore as nuclear density increases (i.e., heavy nuclei or nuclear matter), particle 
masses can drop below the characteristic chiral mass scale and introduce density-dependent 
renormalization in the parameters of the chiral Lagrangian. We have illustrated this with 
BR scaling |14[] but the basic issue can be more general in the sense that it involves relevant 
and irrelevant degrees of freedom as the length scale is varied. This problem needs further 
study. 
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Appendices 

A Renormalization of the VnNN vertex 

Consider the VirNN vertex T^y b (k,q) defined by the kinematics 

N(m N v) -» N{m^v - k - q) + v h {q) + V (fc). 

Restricting ourselves to on-shell nucleons (i.e., v ■ (k + q) = 0) and including the l/m N 
corrections, the vertex has the form 



ttV 



JTT 



^ ^abc^c 

v • q 



9A 



f3 y U 1 f3 

J 77 J 77 



9A 



Sh£(y)-^yl%(k) 

Jir 



2 9 -£s»yh S (y) - ~^q- Shf(y) - 2^y I*(k) - 2^q u I^(k) 
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,9 A 
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2m N / 7r 

r fi,ab 
1 ttV.CT 



(A.l) 



where ^^y CT are the contributions from the counter-term Lagrangian, y = v ■ k and 
h \ (y) = I [^0,4(2/) T ho^(-y)] which were given explicitly in j/J. The functions l\ s (k) = 
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[in(k) =F and I^(fc) are denned by 
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I^ u (k) 



S ■ I (21 + fc)A« 
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d = 4 — 2e is the spacetime dimension and is the renormalization scale. For small 
momentum, they can be simplified to 



I»(k) = -^2Ls l "-(S f "vk + v^S-k)r ] + 0(k 2 ) 
I^(k) = ~^v» S v --e^v a k p r, + 0(k 2 ) 



(A.3) 



where 0(k 2 ) denotes terms equal to, or higher than, second order in fc^ and eoi23 = 1- With 
v ■ k = 0, they become 
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and the divergent quantity r\ is defined in eq,(^5|). 

All the divergences in eq.(A.l) can be removed by the counter-term Lagrangian (p9|), 
that is, the terms with c's, whose contribution to the VirNN vertex is 
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Here cf are finite, //-independent renormalized constants. When v ■ k = v ■ q = we have a 
simpler expression: 

Kv( k ^) = -ta b cT c ^S»-5 ab ^e^v a kpq u \f (k 2 )-fo(0)+c? 



1 9 A 



■ Cg 



Aba qpV pa (k a S*-gfck-S)- A ab (S» k a -g£k- S)V a(i 



Terms proportional to do not figure in our calculation and so are omitted. Note that the 
above equation is finite and //-independent and contains two parameters, Cj and eg which 
were fixed in our calculation by resonance saturation. 

B Renormalization of two-pion-exchange contributions 



The two-body vector currents come from the two-pion exchange graphs and one-pion 
exchange graphs with a four- Fermi contact vertex (see Figure ||): 
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\ (v ■ I) 2 (I 2 - ml) [(I - q{) 2 - ml] [(I + q 2 ) 2 - ml] 

-2i<? hc9 f 2 Y,CA 

Jtt A 



+ (1^2) 



S 1 -(l + f ) T^Tf r{ St 




) r£ i» 


(v ■ I) 2 


(Z-f) 2 -m2 






2 

- ml 



+ (l<-2), 
-i(n x r 2 ) 2 —j 



(u • O 2 (2/ + Q2 ~ QiY 



4/4 7; (Z 2 - m 2 ) [(/ - 91 ) 2 - mj] [(/ + q 2 ) 2 - m 2 } 



J» m (g) oc v-SS" = 



(B.l) 



where 



la 1/3 (I + q)i 



l (v ■ I) 2 (I 2 - ml) [(I + q 2 ) 2 - ml] 



35 



= ^{9auq/3 + gfBuQa- 9a^qu) fo(q 2 ) + Y^K^q 2 )^^ + (v's). (B.2) 

Here and in what follows, (f's) generically stands for terms proportional to v a , vg or v„. 
The function fo(q ) was defined in eq. dAH) and the Ki(g 2 ) was given in @, 



KM 



dz In 



1 - z(l - z) 
-z(l - z)q 2 



mt 



/o m£ - z{± - z)q A 
Since for the Ml transitions, it is sufficient to retain only the terms linear in k, we 



define 



A direct calculation leads to 



8 W 



dk h 



k=0 



(B.4) 



_t 4m ( a )> 



(B.5) 



with 



*&"(<Q = "(2ri + 2r 2 - in x r 2 ) 2 || (s?S?SgS$ - S%S?SgS2) F a0 (q) + (1 ~ 2) (B.6) 
and 

Jl ■ i,- \y~ - hi,- j iyi t Hr - 

rl „2 

(B.7) 



la lf3 

h (v ■ I) 2 (I 2 - ml) [(I + q) 2 - m 2 } 

= -w/o(9 2 ) + |y / dz— _ + ( w 's) 

air* Jo m~ — z[l — z)q z 

where (Vs) denotes terms proportional to v a or vp. 

The divergencies appearing in the one-loop graphs, eq.(B.l), are removed by the 

counter-term contribution coming from the terms with G's and Erf in the counter-term 

Lagrangian eq. (|39[) . Summing up all the two-pion-exchange and counter-term contributions, 

we have 

'A d 



J^(2tt) 



4# dq u 



i[S?, Si-q] f (q 2 ) 



+ (2n + 2r 2 - i Tl x r 2 ) z |! (JL-F^sufr, q) - s?s?s$sZF a/3 ( q fj I 



- (/x <-» v) + (1 <-» 2) 



(B.8) 



where 
^(5i,5 2 ;g) 



S^SiS^Sz — Si Sf S2 Fapvici), 



{g av qp + gpuq a - g a pqv) K (q 2 ) + ^^(g 2 ) ^^ + (Vs), 



16vr 2 



^o(9 2 ) + 



87T 2 Jo W 2 — ^(1 — 2) 



+ (u's) 



(B.9) 



36 



and J^ itc (2ir) is the sum of the divergent part of the loop graphs and the counter-term 
contribution, 



j£L(27r) = G §(n+T 2 y(sr+sn-G%(n-T 2 y(sr-sn 

- \ E% R (n x r 2 y (SrS v 2a - SrS^) . (B.10) 

C Resonance- exchange contributions 

In this appendix, we write down an effective chiral Lagrangian that contains the 
A(1232) and vector mesons [p^ 1 and w M ) in addition to n and N and apply it to the exchange 
currents. Instead of writing a Lagrangian with all low- lying resonances in full generality, 
we shall write only those terms relevant to nuclear processes which can be compared with 
successful phenomenological model Lagrangians employed in nuclear physics. (The axial 
field a\ can be incorporated readily but we will not consider here.) 

In heavy fermion formalism, including the "l/m N " terms, the relevant Lagrangian is 



— f 2 
£ res = B (iv ■ V + 2g A S ■ iA) B + f*Tr{iA„ iA") + ^m 2 Tr(£ + E*) 

1 



2m N 



B ( -V 2 + (v • V) 2 - g 2 A (v ■ iA) 2 - 2ig A {v ■ iA, S • V} ) B 



2m N 



+ m 2 p Tr - iir M ) + l -ml (u^ - * - ^Tr {p^ v ) - \uj^ 

+ C A + C an (C.l) 



with V p = d n - igpPn - \g^^ Puu = d^p u - d v p p - ig p [p p ,p„], 
C A = f n (5m A -iv-V) T M 

and 



V" (-CiA^-i^p^S") B+ h.c. 



(C.2) 



Can = 9^ Va ^ ■ 8 p p u d a UJp + ■■■ (C.3) 

In (|02D , 5m A = m A - m N ~ 293 MeV, TXT = f jH Xl l T mkl and TXB = f jld X^B k e 3lm . 
From the A decay width, we get C = —1.73. The value of has been determined to 
be pt = —7.7 ± 0.5. The anomalous parity term ( |C.3j ) is the same as in [HI] withp^j 

(c ' 4) 

The ellipsis stands for anomalous parity terms that are irrelevant in our calculation. 



* 9 In |Q, Fujiwara et al considered anomalous parity terms in 17(3) with g^ — g p , with g^p-rr = — 
For the ui — > n-y decay, their result is the same as ours. However for p — > n-y, they predict T(p — > n-y) - 
T(u> — > 7T7), which is not in agreement with the experimental data. 
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We determine g p by the KSRF relation, m p = 2g 2 f 2 or g p ~ 5.85, and g w by assuming 
the ideal mixing, g w = 3g p ~ 17.55. To see how good the anomalous parity Lagrangian is, 
consider the decay widths of u — > irj and p — * irj. Our Lagrangian predicts 



a 



r(o> -» 7T 7 ) = — 



9u>pn 



O 



r(p -► vr 7 ) = — 



to . 



to. 



1 



TO, Z 



~ 0.758 MeV, 



TO 



1 1 ~ 0.080 MeV, 



TO, 



(C.5) 



(C.6) 



which are in reasonable agreement the experimental values |35|], r exp (w — > 7r7) = 0.717 ± 
0.050 MeV, r ex P(/9 ± -» 7r±7) = 0.068 ± 0.007 MeV and T ex P(p° -> vr°7) = 0.121 ± 0.031 
MeV. 

The characteristics of this Lagrangian can be summarized as follows: 



It has chiral symmetry and other symmetries consistent with low-energy QCD. 

It respects the vector-meson dominance, universality etc. and satisfies the KSRF 
relation. 

The chiral Lagrangian containing pions and nucleons only are recovered if the m p and 
to w masses are taken to infinity. 



From the Lagrangian (C.l), one sees that only six graphs given in Fig. |7| can make 
non-zero contributions to the vector current, other graphs being either of order 0(m~ 2 ) or 
vanishing. They have simple physical interpretations. Figure 0a is the seagull graph with 
no form factors since ai-meson is absent. Figure |J> is the "pionic graph" with a form factor 
associated with the p- meson dominance. Figures and stem from C an , with Figure 
contributing to the baryonic current. Figure [?]e is new: although it gives a non-vanishing 
short-ranged Sachs moment, it cannot contribute to the n in two-body systems. Finally, 
Figure [7|f is a A-resonance contribution. The explicit forms of the currents derived from 
these graphs are: 



•-4m ( a ) 

^m( e ) 
■&(/) 



-,9 a 



So 



f4 m2 



^ + (i 

"92 



-2i{r\ x r 2 ) 



:9a 



to; 



2), 

Si -qi S 2 - q2 



k^ k • q 



k 2 



rat 



q\ to 2 



to; 



-IT. 



9A 



m. 



92 • S; 



16n 2 f 2 m 2 — k 2 to 2 
QA rn 2 

IT\-T2 



2 9 

qf m{ 
9 2 P 



2 pvaP i. 



9l 



92 • S 2 



Wtt 2 / 2 to 2 — k 2 to 2 — q 2 to 



-i(ti x r 2 ) r-^ 2 



i 

2 TO 



to; 



H-rCgA 



p 

92 • #2 



k 2 



TO 7 



g 2 to 2 



4TO N 0TO A / 2 TO 2 



9 2 2 




+ (1 2), 
^^9/9 + (1^2), 



m 



/9a, 
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p p 



7T 7T 
— — 

P 



(b) 



A 



LO 7T 
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V- — 



(d) 



(/) 



Figure 7: Graphs that contribute to the two-body vector currents. 



+ (l<-2), (C.7) 



where the projection operator 'P"' 3 is defined by 



Vf = v a vP - g al3 - -S?S?. 

o 

We have dropped terms of order 0(m~ 2 ) and order • 0{m~ 1 )^ while keeping terms 
up to 0{m~ 1 ) for the space component and up to 0(m^) for the time component. The 
corresponding magnetic operators then take the form 

g\m n /2 > (x) ~ (x) N 



Kb) 
H(c) 



e T6^ [3 s Tt I yi(X7r) ' 

OA™* f 2 rf,(x) 



16/2 V3 



- 3 yo (x.)] - 4 X) yi(^) ) > 



64n 2 f 2 — \3 

+ + [ m J y 2 (^) - j/ 2 (m,r)] ) , 



/z(d) = - 



eg a 



9 2 P 



32n 2 f 2 m 2 — m 2 \3 



(+) 



mly (x n ) ~ m 3 y (m p r 



+ 4 +) 



mly 2 (x n ) - mly 2 {m u) r) ) , 



//(e) 



o, 



"36m N <5m A /2 \3 ^ s 
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+ (4 +) +4 _) + ^ x) )^)) 



(C.8) 



where yo(x) = yo(x) — 5{x) and 




= Ti • T 2 (<T! + <7 2 ), 



(C.9) 



As emphasized in the main text, the short-ranged contribution /x(e) (i.e., the four-Fermi 
contact term) is equal to zero. 
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